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A lithium dopant in a cuprate spin ladder acts as a vacant (non-magnetic) site accompanied by an 
extra hole bound by a Coulomb potential. We find that, although the undoped ladder spin gap is not 
essentially altered by Li doping, a dopant-magnon bound state appears within the gap. Contrary 
to the case of Zn-doped ladders, we predict that, in the Li-doped ladder, the spin liquid character 
is very robust against antiferromagnetism. We also predict the spatial dependence of the density of 
states in the vicinity of the dopant which could be measured by local spectroscopic probes. 

PACS: 75.10.-b, 75.50.Ee, 71.27.-fa, 75.40.Mg 



Impurity doping of the high-Tc cuprate superconduc- 
tors is an effective tool for exploring the low temperature 
physics of these strongly correlated electronic systems. 
The similarities or differences observed upon doping non- 
magnetic zinc (Zn) and lithium (Li) ions must find their 
explanations both in the strongly correlated nature of 
their host and in the peculiarities of each dopant. For 
example, when doped into the antiferromagnetic (AF) 
phase of La2Cu04, Li [Q is far more effective at sup- 
pressing the Neel temperature than Zn ^ , although both 
enter the same planar Cu(2) site. However, Li, by releas- 
ing an extra hole (due to the difference of formal valence 
of Cu^"*" and Li"*") can theoretically lead to an essential 
perturbation ||^ while Zn was shown to only slightly en- 
hance the AF correlations in its vicinity Q. 

Strikingly, despite different behavior in AF systems, 
non-magnetic Zn and Li ions behave quite similarly in 
conducting YBa2Cu306+a; (YBCO) by inducing local 
magnetic moments |^]. This suggests that, in the case of 
Li, the extra hole is not bound, unlike La2Cui_a;Li2;04. 
The induced moment sits predominantly on the four 
nearest neighbor (NN) Cu of the dopant and exhibits 
static lH] and dynamic susceptibilities reminiscent of 
a Kondo-like behavior with a large range of Kondo tem- 
peratures. In fact, screening of the induced moment by 
a conduction hole resulting in an impurity-hole bound 
state was predicted theoretically |^ using a description 
in terms of a vacant site embedded in a correlated t-J 
host. In superconducting Li-doped YBCO the moment 
was found to be weakly modified below Tc Possible 
reduction of screening due to pairing was also investi- 
gated theoretically pO| ]. 

Spin ladders [0 offer a perfect system to investi- 
gate impurity doping in a spin liquid or Resonating Va- 
lence Bond [|l2| (RVB) state, i.e. a state with short AF 
correlations and a spin gap. Clearly, a reliable micro- 
scopic description of dopants in ladders would greatly 
contribute to the understanding of their two dimensional 
(2D) analogs as well. In fact, Zn doping into the spin- 
1/2 Heisenberg two-leg ladder SrCu203 leads to local 



moments (as in metallic YBCO) and to an instability 
of the gapped spin-liquid state towards an AF ordered 
state |13| at low enough temperature. Formation of local 
moments as well as a rapid suppression of the spin gap 
was obtained theoretically in a Heisenberg ladder using 
the vacant-site impurity model [Q. Further numerical 
simulations have led to an effective model of local 
moments interacting via an exchange interaction which 
rapidly decays with the impurity distance. 




FIG. 1. Schematic representations of a Li-doped spin lad- 
der. The cross (resp. the circle) stands for the impurity 
site (resp. injected hole) and the thick lines corresponds 
to the attractive potential, (a) Singlet GS where spins are 
paired up into spin singlets (shown as shaded areas). Sites 
are labelled for convenience, (b) Lowest triplet excitation: 
dopant-magnon bound state. 

Doping Li into a two-leg spin ladder is still an open 
and fascinating problem both experimentally and theo- 
retically. Novel physics can be expected due to the ad- 
ditional (w.r.t the case of Zn) intrinsic doping associated 
with the substitution of a divalent (Cu^^) by a mono- 
valent (Li+) ion. In this Letter, we address this issue 
using a specific vacant-site model for Li+ embedded into 
a correlated host. We show that Li doping does not in- 
troduce low-energy spin excitations as in the case of Zn 
but instead leads to the formation of a magnon-impurity 
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bound state (BS) just below the spin gap of the undoped 
ladder. This suggests that, unlike Zn-doped ladders, Li- 
doped ladders would keep a robust spin liquid character 
at low temperature. 

For dilute concentrations, a single dopant in the geom- 
etry of Figs. |l| suffices. We model the substitution of a 
Cu^"'' ion by a Li+ atom by a vacant inert site which is 
a hard-core potential for holes. Because of the different 
valences of Cu^+ and Li+ a hole is injected into the spin 
ladder and it feels a static attractive Coulomb potential 
centered on the Li+ site. A realistic Hamiltonian reads: 

= J ^ (S, . S, - in,n,) (1) 

<ij> 

<i,j>.,cr I] 

where the notations are standard and the prime means 
that the sum over NN links < ij > is restricted to the 
bonds not connected to the dopant site. We have mod- 
eled the attractive hole potential by a static NN potential 
—V as depicted on Figs. 0. The summation over Ij only 
runs over the NN sites of the impurity. Note that, for 
simplicity, we restrict ourselves to the case of a magnet- 
ically isotropic ladder i.e. with equal rung and leg cou- 
plings, J. Numerical results discussed in the following are 
obtained by Exact Diagonalisations (ED) of small peri- 
odic ladder rings (up to 2 x 12) supplemented by Density 
Matrix Renormalisation Group (DMRG) calculations on 
larger open systems (up to 2 x 64). In open ladders, the 
impurity site is placed in the center. 
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FIG. 2. Various binding energies to the impurity site (see 
text for definitions) vs J/t obtained by ED on a 2x11 lad- 
der. The amplitudes V used for the attractive potential are 
indicated on the plot, (a) Binding energy of the single hole 
in the singlet GS. For comparison, half of the hole-pair bind- 
ing energy is also shown, (b) Binding energy of a hole to the 
impurity-magnon system (see [Q) in the lowest triplet exci- 
tation, (c) Binding energy of a magnon to the impurity-hole 
system (see [Q) in the lowest triplet excitation. 



First, we investigate the localization of the injected 
hole vs the strength of the Coulomb potential. Following 
Ref. we calculate the so-called hole-dopant binding 
energy defined as, 

Af=°.i, = Eo{lh,U)+EoiOh,Oi) (2) 
-iE„{lh,Oi) + EoiOh,U)), 

where Eq^u h, m i) is the GS energy of a ladder with n = 
or 1 (m = or 1) holes (dopants). This quantity is nega- 
tive when the hole and Li ion form a stable bound state. 
Since Li-doping results in the simultaneous removal of 
two spins we expect a magnetically inert state i.e. a sin- 
glet {S — 0). As seen on Fig. ||(a) a stable BS is found for 
almost all couplings, even in the case = 0. It should 
be noticed that the magnitude of the binding energy is 
slightly larger than in the case of a 2D planar geome- 
try ^ . Fig. |^(a) offers a pictorial representation of such 
a state where the absence of any local moment can be 
clearly understood from the RVB nature of the host (all 
remaining spins are paired up in singlets). Although, 
our calculation shows that a single Li-dopant binds the 
injected hole in its vicinity, caution is required at a finite 
concentration. In this case the possibility of other "de- 
cay channels" has to be considered, e.g. 2 holes from 2 
dopants recombining into an itinerant hole pair p6[ | . This 
possibility is ruled out by the fact that the dopant-hole 
binding energy is always smaller than half of the hole-pair 
binding energy Aahoioa calculated on an undoped ladder 
as seen on Fig. |^(a), even when V ^ 0. At low con- 
centrations we then expect decoupled quenched dopants 
each binding one hole. Since the spatial extensions of the 
isolated BS's are quite small (typically ^ = 2 to 4 rungs 
even when V — 0), the system is expected to remain 
insulating up to large doping. 

We now turn to the magnetic properties of the Li- 
doped spin ladder in order to compare with Zn doping. 
For this purpose, we compute the dynamical spin struc- 
ture factor 5(q, ui) as measured in Inelastic Neutron Scat- 
tering (INS) experiments for an undoped ladder (shown 
on Fig. ^(a)), a Zn-doped ladder (modeling Zn as a neu- 
tral vacant site and shown on Fig. ^(b)) and a Li-doped 
ladder described by Eq. (|^) (shown on Fig. ||(c)). Note 
that S{c[,uj) quantifies the spectral weight given by the 
matrix element \{n\S^\0')\'^ between the singlet GS |0) 
and the magnetic excitation |n) at energy co. The RVB 
picture, in which spins are paired up into short range 
singlets (for example on the rungs as in Figs. ^), turns 
out to give a nice qualitative understanding of the mag- 
netic properties. In the undoped ladder, a triplet excita- 
tion (magnon) is well described by exciting a rung singlet 
into a triplet. Fig. ^(a) clearly shows the single magnon 
dispersion with a minimum at q = (tt, tt) which defines 
the spin gap ijl^ of the undoped ladder. Introduc- 
ing a Zn atom on a rung releases a free spin- 1/2 which 
leads to zero-energy spin fiuctuations predominantly at 
q = (tt, tt) as seen in Fig. ||(b) while the undoped ladder 
magnon mode still survives. When two Zn dopants are 
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FIG. 3. Spin structure factors 5(q,tj) calculated on 2 x 10 
ladders. Data for q± = n only are shown. From bottom 
to top, the different curves correspond to decreasing gy mo- 
menta, from q = (7r,7r) (bottom) to q = (0, tt) (top). For 
clarity, reducing scaling factors (as indicated) are applied on 
some curves, (a) Undoped periodic ladder; (b) Spin ladder 
doped with a single Zn impurity (full line) or two Zn im- 
purities separated by the maximum distance on the same leg 
(dashed line) ; (c) Ladder doped with a single Li impurity with 
t = 2J and V/t = 0.5. 



introduced, our data, in agreement with Ref. |14|, show 
that the two local S=l/2 moments experience a weak ef- 
fective exchange interaction, Jgfj, which decays rapidly 
with separation. In that case, a small spectral weight at 
q = (tt, tt) and low energy uj ~ Jgfj below the undoped 
spin gap is seen. The case of Li-doping shown on Fig. ||(c) 
is drastically different with no weight at small energy. 
Indeed, since a Li+ dopant introduces simultaneously a 
bound hole, it does not release a free spin as for Zn but 
rather produces a new type of excitation located just be- 
low the unperturbed spin gap, namely a collective bound 
state of a (let us say rung) magnon with the hole and the 
Li ion as naively depicted on Fig. ^b). Its binding energy 
defined as the energy separation w.r.t. the free magnon 
energy Ag remains in general quite small (in absolute 
value) as seen directly in Fig. |^(b) (and qualitatively in 
Fig. 11(c)). Therefore a drastic reduction of the spin gap 
does not occur in this case. 

The physical origin of this BS is of particular inter- 
est. Simple qualitative arguments as well as extensive 
numerical data show that the origin of binding in this 
channel can be attributed to the hole kinetic energy gain 
associated with the spin polarization. DMRG calcula- 
tions show that the spatial extents of the charge and 
spin perturbations associated with the BS are quite dif- 
ferent. While the hole is localized on a scale of a few 
rungs from the dopant as seen in Fig. ^a), the distri- 
bution of the rung magnetization can extend to large 
distances as seen in Fig. ^(b). This is in agreement with 
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FIG. 4. Hole rung density (a) and Sz rung density (b) 
along the ladder direction in the lowest energy triplet state 
calculated by DMRG for J/t = 0.5. The rung density is 
defined as the algebraic sum of the densities (if any) on the 
two sites of a given rung. Different values of V/t (as indicated) 
are shown. 



our other finding that the energy costs (i.e. the absolute 
value of the related binding energies |lq]) of the two vir- 
tual decay processes leading either to a free hole and a 
neutral dopant-magnon complex or to a free magnon and 
a dopant-hole complex are quite different in magnitude 
as seen in Figs. ||(b) and ^(c). In fact, the attraction 
—V has the opposite effect on the two quantities; while 
increasing V binds the hole more strongly, it also lim- 
its the ability of the hole to reduce its kinetic energy by 
moving in the spin polarized background of the magnon, 
and so the binding energy of the magnon to the dopant- 
hole complex is reduced |l^]. From the analysis of the 
behavior of the magnetic size of the BS (deduced from 
Fig. ^(b)) or alternatively directly from its binding en- 
ergy plotted on Fig. ^(c), we conclude that above a criti- 
cal value of V (typically one gets Vc/i — 2 for J/t = 0.5) 
the magnon escapes from the dopant. 

For completeness, we calculate also the local density of 
state (DOS) in the vicinity of the dopant. Since this DOS 
might be accessible experimentally by Scanning Tunnel- 
ing Spectroscopy it can test the validity of the model. 
Results are shown in Fig. ^ corresponding to spatially 
resolved DOS spectra in the vicinity of the dopant. The 
UJ > {uj < 0) spectra give the weights of the neutral 
(charged) target 5' = 1/2 states which are accessible by 
removing the hole (adding an extra hole) in the singlet 
GS. The large peak seen at vanishing positive energy on 
site number 1 i.e. on the same rung as the dopant corre- 
sponds to a resonant magnetic spinon-dopant BS. Other 
spinon-dopant resonances of higher energies are also seen 
with larger spatial ranges i.e. spread over two or three 
sites of the leg opposite to the dopant. Similarly, w < 
resonances are seen when adding an extra hole on the 
same leg as the dopant e.g on sites number 3 and 5. The 
lowest resonance energy (in absolute value) is obtained 
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FIG. 5. Local DOS around a Li dopant obtained by ED 
of a 2 X 9 ladder with J/t — 0.5 and V/t — 0.5. Each panel 
corresponds to a site in the vicinity of the dopant (site labels 
correspond to those of Fig. |l|(a)). Occupied (empty) elec- 
tronic states are shaded (left blank). Note that, for a dopant 
concentration x, the total w > (resp. cj < 0) integrated 
weight scales as x (resp. | — f ). 



when the second hole is added on sites number 3 next to 
the dopant. Note also that the local a; > (w < 0) inte- 
grated weight provides directly the local hole (electron) 
density in the GS. Hence our data show that the bound 
hole is located mainly on the opposite leg to the dopant 
and extends roughly over three sites. 

Our theory can be directly tested, if Li can be sub- 
stituted for Cu in the ladder compound, SrCu203. The 
extra bound hole around a Li dopant should ensure that a 
free local S=l/2 moment is not created in contrast to the 
case of Zn doping. Hence we predict that SrCu2-a;Lia;03 
should not order antiferromagnetically at low tempera- 
ture unlike SrCu2-a;Znj:03 (ReL |13|). Further the na- 
ture of the magnon-dopant BS, the charge distribution 
and local DOS could also be examined experimentally 
in these systems. However, our analysis raises interest- 
ing questions regarding the close similarity between Li 
and Zn substitution reported in superconducting YBCO 
samples. In particular if we interpret the spin gap phase 
in underdoped YBa2Cu306.6 as a doped d-RVB phase, 
then there should be a close similarity to the behavior 
of the doped ladder. However Bobroff et al. ^,|| re- 
port a free S=l/2 appears for both Zn and Li doping of 
these samples. A possible way to reconcile this apparent 
contradiction is to postulate that Li+ does not bind a 
hole in YBa2Cu3_a;Li2:06.6j but does in La2Cui_a;Lia;04. 
This could occur if the mobile 0^^-ions in the chains 
were repelled from the Li^-ions on the planes. A test 
of this hypothesis can be made by doping Li"*" and Zn^"*" 
in YBa2Cu408 which as a stoichiometric compound has 
no mobile 0-ions. Our analysis then predicts free S=l/2 
moments only for Zn-doping and not for Li-doping in this 
case. 



In conclusion, our analysis predicts a clear distinction 
between the magnetic properties of the two non-magnetic 
ions, Zn^+ and Li+ when doped into spin liquids due to 
the binding of a hole in the latter case. Experiments to 
test these predictions are proposed. 
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